Stochastic Representation of Deterministic Interactions 
and Brownian Motion 
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Exact generalized stochastic representation of deterministic interaction between two dynamical 
(quantum or classical) systems is derived which helps when considering one of them to replace 
another by equivalent commutative (c-number valued) random sources. The method is applied to 
classical Brownian motion of a particle in a gas, and statistics of this motion is reduced to statistics 
of the gas response to perturbations. 
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I. INTRODUCTION 

In recent years the old variety of problems about in- 
teraction between one or another "dynamical system" 
(DS) and its environment (named also "thermal bath" 
or "thermostat" ) have attracted new attention as a part 
of the quantum computation problem. Frequently, the 
environment contains too many routine degrees of free- 
dom to be of special interest in itself, and one would 
be satisfied by a simplified statistical description of DS 
in terms of its own degrees of freedom only but with 
adding of random sources, or "Langevin forces", which 
effectively substitute the "thermostat". In this respect, 
during last century many different approaches were pro- 
posed. Perhaps, the two extreme points among them are 
the highly time-nonlocal description of the DS evolution 
by means of path integrals and influence functionals 0, 
on one hand, and the mathematical theory of quantum 
semigroup evolution operators corresponding to purely 
time-local kinetics 0, on opposite hand. 

In our previous works [3- Q the new intermediate 
approach was suggested, so-called "stochastic represen- 
tation of quantum interactions" . Its peculiarity is that it 
involves time-local stochastic evolution equation for den- 
sity matrix of DS, with commutative (formally c-number 
valued) random sources which generally undergo time- 
nonlocal correlations. In contrary to phenomenological 
"Langevin equations" for DS variables, this equation is 
exact, in the sense that quite definite and transparent re- 
lations are established between complete statistics of the 
random sources and private statistical characteristics of 
the thermostat perturbed by classical (c-number valued) 
forces. These relations ensure adequate construction of 
approximate models of the thermostat noise in its con- 
nection with the thermostat induced dissipation. 

In the present paper, most general formulation of 
this approach is presented which unifies both quantum 
and classical mechanics, at both Hamiltonian and non- 
Hamiltonian interactions between DS and thermostat. 



The results are illustrated by Brownian motion of clas- 
sical particle in classical gas, the "trivial" phenomemon 
which still remains far from complete consideration (as 
well as Brownian motion of gas molecules themselves, see 
and references therein). 



II. STOCHASTIC REPRESENTATION OF 
TWO-SIDE INTERACTION 

Let p{t) be joint statistical operator (distribution func- 
tion, density matrix) of combined system "DS plus ther- 
mostat (environment)" and L(t) its evolution operator, 
so that 

dp{t)/dt = L{t)p{t) (I) 
We suppose that L(t) has the bilinear form as follows: 



Lit) — Ld(t) + Lb + L int , L int — ^ A^A^ 



(2) 



where operators Ld(t) and Lb are responsible for au- 
tonomous evolutions of DS and thermostat ("bath"), re- 
spectively, while the pairs and A* for their interac- 
tion. The indexes "d" or "b" mark objects relating to 
DS or thermostat, respectively. Notice that all of Ld and 
A^ commute with all of Lb and A„, since these two sets 
of operators act onto variables of two different systems. 

In general, Ld{t) contains a part which describes ex- 
ternal observation of DS. For instance, if DS undergoes 
quantum Hamiltonian evolution, then 



Ld(t)p = v k{t)Jk °P + i[p, H d (t)]/h 



(3) 
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where Hd{t) is Hamiltonian of DS; Jk are some opera- 
tors belonging to DS and representing quantum variables 
under observation; Ufc(t) are test functions related to Jk ; 
[, ] and o mean commutator and Jordan symmetric prod- 
uct, respectively: 

[A, B] = AB — BA , Ao B = (AB + BA)/2 
In case of classical Hamiltonian evolution, similarly, 

Ld{t)p = Mt)JkP + {p, H d (t)} , (4) 
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with {,} being the Poisson bracket, {A, B} = f^^f- 

W T§> (l ano " P are canon i ca l coordinates and momenta 
of DS), and Jk some phase functions (i.e. functions of 
q and p ) under observation. Then the trace, Tr p (or, 
in classical variant, integral over total phase space), dif- 
fers from unit and gives characteristic functional of the 
observables Jk(t), that is generating functional for their 
multi-time correlators [HEllil]- 

Remembering that L d and commute with L b and 
A* , it is not hard to justify the formal identity 



cxp 



L d (t)+L b + J2<K 



dt 



(5) 



3Xp 



L d (t) + Va„(t)A^ dt\ | a(t)=0 



2Xp 



where exp denotes chronologically ordered exponent. 
Because of the chronological ordering, this formula is 
the mere consequence from the identity exp(Cq02) = 
exp(Oi^) exp(aC?2) |a=o > with 0\ and O2 being two 
mutually commuting operators ([OijCy = 0). 

From the other hand, if a n (t) is a set of random pro- 
cesses with a given characteristic functional 



*{/} 



exp 



^ fn{t)a n {t)dt 



(the angle brackets mean statistical averaging), then the 
average of arbitrary functional ${a} can be formally rep- 
resented as 



{ — \ | Q=0 



(6) 



Now, let us assume that more or less far in the past 
the joint statistical operator p(t) was factorized: p(to) = 
Pdo*pbo (for instance, at to = —00). Then comparison 
between the Eqs. |5] and clearly results in two con- 
clusions. First, the partial statistical operator of DS, 
~p d (t) = Tr b p(t) , can be represented in the form 



p d (t)=Tx b p(t) = (p d (t)) , 



(7) 



where Pd(t) is a solution (under initial condition p d (to) 
= pdo) to the stochastic evolution equation 

dp d {t)/dt = ^L d {t) + J2 a n (t)K}pd(t) , (8) 

with a n (t) being random sources, and the angle brackets 
() standing for the average with respect to these sources. 
Second, characteristic functional of these sources is de- 
fined by 



Tr 



b cxp 



L b + Y,fn(t)A b n 



dt 



PbO 



Hence, from the point of view of DS its interaction with 
thermostat is equivalent to its separate evolution but dis- 
turbed by random sources and therefore described by 
stochastic statistical operator of DS, p d (t). The aver- 
age (pd(t)) , being a functional of the test functions Vk(t) 
introduced in J3J and exactly reproduces not only 
current statistictical state of DS but also all the multi- 
time correlators of the observables Jk- 

The Eas l7l9l give what we call stochastic representa- 
tion of interaction between two systems. In this repre- 
sentation, the DS and thermostat evolutions can be anal- 
ysed separately. The cost of such possibility is that both 
DS and thermostat should be considered under arbitrary 
time-dependent perturbations (a n (t) and f n {t), respec- 
tively). The examples 0, 0, demonstrate that this 
cost may be not too large, since when considering DS we 
become saved of many degrees of freedom of thermostat 
and instead must deal with a few random sources only. 
Importantly, like c-numbers, the sources always commute 
one with another and with any other objects. Besides, 
possibly, their exact statistics which follows from @, can 
be replaced by some simplified (semi-phenomenological) 
statistical model. Dependently on concrete contents of 
operators A^ and A n , the sources a n (t) behave either 
as literally classical random processes or as commuta- 
tive "ghost fields" possessing unusual statistical proper- 
ties (see 0, 0, and below) . 



III. BILINEAR HAMILTONIAN 
INTERACTIONS 

Let a quantum system "DS plus thermostat" has the 
Hamiltonian 

H(t) = H d (t) +H b + H rnt , H int =J2 B 1* D i ' ( 10 ) 



whose interaction part, Hi nt , is bilinear, with operators 
Dj and Bj acting in different Gilbert spaces of DS and 
thermostat, respectively, and therefore commuting one 
with another. The corresponding joint evolution operator 
is L(t) = L d (t)+ L b + L int , where L d (t) is given by ®, 
L b p = i[p, H b ]/h , and 



(11) 



Because of [Dj,Bj] = 0, for any pair Dj and Bj the 
equalities take place: 



[p, D * B] = [p, D] o B + D o [p, B] 



[p oB,D] + [po D, B] 



(12) 



exp 



^fn{t)a n {t)dt 



(9) 



Hence, the interaction part of the evolution operator has 
just the bilinear structure as in J2J). We see that every 
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term Dj*Bj from (|10|) produces two terms in Li nt in J3J, 
namely, 



Lint - A 1i<? A % , 

j a=l,2 



(13) 



with the terms defined by 



A*0 = i[0,D 3 }/h, A'Oe^oO, (14) 



K d j2 = DjoO, A b j2 = i[0, Bj]/h 

(here O is an arbitrary operator). 

Correspondingly, every term of Hi nt results in two ran- 
dom sources in (JSJ), aj a (t) [a = 1,2). It is convenient 
to rename them as ctji(t) = Xj(t) , ctj2(t) = Uj(t). After 
that the Eq|Hl takes the form 

dp, 



For details and another equivalent representations of the 
characteristic functional see U3- LH- 0] It should be empha- 
sized that due to the causality principle any of the cor- 
relators l|18l) turns into zero if max m {r m } > maxj jij}. 
In addition, without loss of generality, the operators Bj 
can be defined in such a way that (xj(t)) = . 

In case of thermodynamically equilibrium thermostat, 
the important relations between pair correlators take 
place, namely, 



(Xj(T)x m (0)) = 



cos(ujT)S jm (u)— , 



(x j (r)y m (0))=d(r) x 



' f°° ■ , ^ , a , 

- J sm(wr)tanh ^— j bj m (w) — 



(19) 



_ L d (t)p d + Y Vi (t) D j ° Pd + Y x i (*) ^ D i\ Here T is the thermostat temperature, i?(r) is Heavyside 



According to and (|14|l , characteristic functional of the 
sources in this equation looks as 



(15) function, and Sj m (uj) is a non-negatively defined spectral 



bPb , (16) 



where operator pb (effective partial statistical operator 
of thermostat) is defined as solution to the equation 



dpb T 

Ht =LbPb ~ 



matrix. The connection between the two types of correla- 
tors in (|19|l merely expresses usual fluctuation-dissipation 
theorem. In accordance with l|19|) . if we neglected the 
sources y(t), it would be equivalent to infinite tempera- 
ture. 



IV. NON-HAMILTONIAN INTERACTIONS 

We can obtain serious generalization of the represen- 
tation (0-10 if replace (JT2J! by 



with initial condition pb(ta) = pbo- This equation de- 
scribes separate evolution of thermostat, under its per- 
turbation by classical (c- number valued) forces fj (£) con- 
jugated with the observables Bj, and besides under mea- 
surement of the same observables. The gj[t) are the 
probe (test) functions of the measurement. 

In the stochastic "Liouville-Langevin equation" (|15fl . 
the sources x(t) play the role of usual (realistic) random 
forces (potentials), while y(t) are "ghost" ones. Indeed, if 
we put on g{t) =0 in 117|) then the evolution described 
by these equations becomes purely unitary. Therefore 
Tr b pb = 1, and the Eg ITo! yields (exp / f(t)y(t)dt) = 1. 
Consequently, all the statistical moments of y(t) and 
their correlations between themselves are zeros. Never- 
theless, the joint statistical moments and cross correla- 
tions of y(t) and x(t) are nonzero, being responsible, in 
particular, for dissipation in DS 0- [J Q According to 
IJEJ and 03, 



(20) 



A%0 = D' 3 oO. 



thus involving not two but four different observables at 
any j, with D' and D belonging to DS while B' and B 
to thermostat. Instead of H15|) and l|17fl . we come to the 
representation as follows (summation over j is taken in 
mind) : 



= L d (t) Pd + y(t)D' o p d + x{t) l -[p d , D] , 



dpb T 

dT = LbPb 



g(t)B'op b 



f{t)-[p b ,B] 



(21) 



(22) 



n 



Sf(T m ) 

' f=0 

(18) 

where B(t, /) are the thermostat observables B consid- 
ered as functions of time and functionals of the Hamil- 
tonian perturbation characterized by the forces f(t) . 



This is accompanied by the same definition l|16|l of the 
characteristic functional. 

Obviously, the interaction defined by (|13|l and (j20|l is 
not reducable to a bilinear Hamiltonian interaction like 
l[TU)l. since any pair of the sources (x(t) and y{t)) replaces 
not one but two different observables of the thermostat 
(B' and B). As in l|21|l . y(t) are "ghost" noise sources. 
The formula (|18|) extends to this non-Hamiltonian case 
under substituting variables B'(t, f) in place of variables 
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B(t,f) in p8[). An exact contents of the correlators 
(II B'fyjt /)) i s directly implied by the definition (Jfi) . 

The most general form of non-Hamiltonian interaction 
similar to (|21() and i|22|) is expressed by the formulas 

L in4 = ^{A^' + AfA)} , (23) 



= + $>i(*) A f + Xj (t)Af}p d , (24) 



F(rj), and ^ transforms into ^ — Second, 
introduce the operator 

L b = L s -^F(r,). A (28) 

Formally, L b describes the gas in presence of immovable 
scatterer fixed at zero point. After that i|26|) transforms 
into 

L = L d (t) + L b + L mt , (29) 



dfb 
dt 



(25) 



again accompanied by (|16|l . where A^ and Aj are some 
generators of unitary (phase volume preserving) individ- 
ual evolutions of DS and thermostat, respectively, while 
Aj' and A b - generate non-unitary phase volume exchange 
between DS 



and thermostat 
one for another). 



(their mutual observation 



L i 



dP 



drj 



(30) 



Hence, Li nt has just the structure (|23|l . and we can 
identify the A's (replacing their indexes by natural vec- 
torial notations): 



A d = - 



_d_ 
dP ' 



V. BROWNIAN MOTION 

Examples of the bilinear Hamiltonian interaction can 
be found in 0, 0, [j| . As concrete example of the previ- 
ous non-Hamiltonian case, consider Brownian motion of 
classical particle in classical gas. 



M 



E 



<9r%- 



Correspondingly, the Eqs. 1241 and 1231 look as 



dpd 
dt 



L d {t)+y(t)- V-x(t)- A tf)/l . 



(31) 



(32) 



A. Brownian particle in a gas 

Let R, V, P and M be position vector, velocity, 
momentum and mass, respectively, of Brownian parti- 
cle (BP), t\,, Vj , pj and m be analogous quantities of 
gas particles, and the latters interact with BP by mean 
of potentials U(rj — R). The joint evolution operator 
is merely the Liouville operator of the system "DS and 
thermostat (gas)": 

L = L d {t)+L g + Y J F{r j -R)-{^p-^), (26) 

where F{r) = — VU(r), and L g is the Liouville operator 
of the gas itself. Generally L d includes three terms: 

L d {t) = ik(t) ■ V - V ■ A - Fext (t) . , (27) 

with F ext {t) being an external force applied to BP, and 
k{t) the test (probe) function for observing BP's velocity 
(ik(t) plays the same role as v(i)'s in (0J). 

Let us reduce expression l|26(l to the form of non- 
Hamiltonian bilinear interaction J2J, (|23|l . First, make 
the non-canonic change of variables, namely, consider po- 
sitions of gas particles in terms of the relative distancies 
rj —R and redenote the latters as r,-. Under this change, 
L g preserves its form, but F(rj — R) transforms into 



(33) 

to be accompanied by (|27|) and l|28|) . For the character- 
istic functional l|16|l let us introduce the designation 

"{<?, /} EE (Jl9(t)Mt) + f(tyy(Wt^ = Tr b pb (34) 

Clearly, in (|32|l x(t) is actual random Langevin force pro- 
duced by the gas. In l|33|l . factual physical dimensionality 
of the "force" f(t) is velocity. 

Eventually we are interested in the characteristic func- 
tional of the velocity and displacement of BP, that is 

e^F^jEE^J^^) (35) 

It can be found in two steps: 

e{ik : F ext7 x,y} ee J Pd dRdV , (36) 

0{ik, F ext } = {0{ik, F ext , x, y}) (37) 

The first step is very easy, due to that the coefficients 
of the stochastic Liouville operator on right-hand side of 
Eg 1321 (first-order differential operator) do no depend on 
R and are linear functions of V. Supposing, without loss 
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of generality, that initially BP was fixed at zero point 
at its phase space, R(t ) = 0, V(t ) = 0, i.e. p d0 — 
S(R)6(V) , we can obtain 



(38) 

M 



Q{ik,F ext ,x,y} = 

exp j f [ifc(ta) + y(t 2 )\ ■ [F ext {tx) + x(ti)\ 

In opposite, the second step, that is averaging H38[l with 
respect to random trajectories x(t) and y(t), is rather 
nontrivial problem. 

B. Path integral formulation 

Of course, we can rearrange $ and \& in the identity © 
and besides replace a n by c n a n with c being any nonzero 
coefficients. Therefore let us write 



&{ik,F ext } = 



s s 

(-) { ik, F ext , — , — \ E{ig, /} 



-lg=/=o 
(39) 

Due to simple specific structure (|38|l of the functional 6, 
we can easy transform the latter differential expression 
into integral one if use functional analogue of the formal 
identities 



exp it 



d 2 



dadb 



6(a)6(b) 



1 



2.1TT \ T 



iab 



exp it 



<:) 2 



$(a,b) 



1 

2ttt 



dadb 
da / db exp 



(40) 
(41) 



a=b=0 



The final result of subsequent natural algebraic manipu- 
lations is the path integral 



e{ik,F ext }= / DgDfE{ig,f}x 



(42) 



exp < i 



Mf(t) 



dg{t) 

dt 



k(t)f(t)+g(t)F ext (t) 



dt 



with DgDf = Y[ t [Mdg(t)df(t)/2n] being path differen- 
tial in some reasonable space of g(t),f(t) trajectories. For 
brevity, here and below we omit the central dot symbol- 
izing scalar product. 



C. Limit of ideal gas 



with and N being total volume and number of atoms of 
gas, respectively. Naturally, the initial one-particle dis- 
tribution function can be supposed homogeneous (space- 
independent at least far from BP), for instance, Wo(r, v) 
= Wo(v), while Wo(v) arbitrary even (spherically sym- 
metric) velocity distribution. In particular, the reason- 
able choice for Wq(v) is usual Maxwellian distribution 
corresponding to thermally equilibrium gas with a defi- 
nite temperature. 

In the "thermodynamic limit" N — > oo, f2 — > oo, 
N/Cl — » v (y =const), statistics of the thermostat noise 
gets the general Poissonian form (i.e. becomes infinitely 
divisible) : 



= exp(i/Q{g,/}) 



(43) 



Here functional Q{g, f} characterizes contribution to the 
overall noise from one atom: 

Q{g, f}= J [W(t, r, u) - W (r, u)] drdu , 

with function W(t, r, v) undergoing the equation 
dW { d Fir) d 1 

(44) 

and the initial condition W(to,r,v) — Wo(r, v). Hence, 
analysis of Eqs. 1331 and l3~4l in company with Eq[2Hl re- 
duces to the one-particle problem. 

In EqElthe first right-hand term says that an atom 
is observed by BP. The next two terms say that it is 
dynamically perturbed by BP, as if an atom would have 
the Hamiltonian H = [p — mf(t)] /2m +U(r). Hence, 
the "force" f(t) representing velocity of BP acts like an 
effective time- varying vector potential. 

Alternatively, we can consider this one-particle prob- 
lem in terms of the factual velocity of an atom, u = 
p/m — f(t) , instead of its momentum p or v = p/m. 
Then Eq^| transforms into 



dW 
~dT 



= \-g{t)F{r)-u 



0_ 

Or 



dm F{r) 

dt m 



d_ 
du 



W 

(45) 

In agreement with the above mentioned condition V(to) 
= , we can assume that /(to) = and that initial con- 
dition for the W treated as a function of the factual ve- 
locity u = p/m— f(t) is W(to, r, u) = Wq{t, it). Because 
integration over u in 143(1 is indifferent to shift of the in- 
tegration variable, u — > u + f(t), the function W(t, r, u) 
in the characteristic functional (|43|l can be mentioned as 
solution to EallSI 



In the simplest case, let the gas be ideal, i.e. there are 
no interactions between its atoms and besides no initial 
statistical correlations between them: 



L„ 



d 



N 

fc=i 3 



N 



, Pbo = Yl 



k=l 



W (r k ,v k ) 



VI. DISCUSSION AND RESUME 

The above example once again demonstrates (see also 
0- El 0] ) how the stochastic representation formalism can 
help in reducing many-particle problems to one- or few- 
particle problems. It will be interesting to find out how 
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far this approach conducts in analyzing factual statistics 
of Brownian motion beyond conventional Gaussian ap- 
proximation. 

Of special interest is one-dimensional motion which 
arises in the problem about massive piston sliding in (in- 
finitely long) tube filled with an ideal gas It was 
rigorously proved (see [6j for the review and original re- 
sults) that in the limit of infinitely heavy BP (piston) 
its velocity tends expectedly to the Ornstein-Uhlenbeck 
random process possessing Gaussian statistics. 

From the other hand, our (mathematically non- 
rigorous) attempt Q (see also Q) to reconsider deriva- 
tion of Boltzmann kinetics of classical weakly non-ideal 
gas, in the Boltzmann-Grad limit, demonstrated that 
statistics of Brownian motion (self-diffusion) of gas atoms 
is essentially non-Gaussian because of low-frequency fluc- 
tuations in diffusivities of the atoms (uncertainty of the 
diffusivities). 

Hence, the task naturally arises about non-Gaussian 
corrections to statistics of piston's motion at large but 
finite ratio M/m . In principle, it is sufficient to as- 
sume symmetrical situation when initially both "left- 
hand" and "right-hand" gases are in identical (mirror 
symmetrical) states. In order to keep constant piston's 
velocity relaxation time, the gas density v (number of 
atoms per unit tube length) should be maintained pro- 
portionally to M/m . 

Due to the symmetry, in Eg 1431 

Q{gJ} = Q+{gJ} + Q+{-g,-f} , 

where Q+{g, /} is connected to solution of Eq^Jor Eg 1431 
for the right-hand gas atom by the same Eq031 The 
source of possibly non-trivial corrections (as well as non- 
trivial mathematical difficulties, see 0) is the recollisions 
between piston and atoms. 

Indeed, in an infinitely long container (tube), let before 
a collision occuring at t = the piston's velocity, V(t) , 
and a "right-hand" atom's velocity, v(t) , were equal to 
Vb and vq , respectively, of course, satisfying vq < Vq . 
After collision the atom's velocity equals to v = 2Vq — 
vq, while the distance, 5R(t), between piston and atom 



changes as 

SR(t) =vt- R{t) , R(t) = [ V{t')dt' 

Jo 

The recollision will happen if 5R(t) — at t > 0. Sup- 
pose that after the collision the atom runs on the left, i.e. 
v = 2Vq — vq < , which means that it was overtaking 
the piston, i.e. Vq < . Besides let V(t) is symmetric 
Ornstein-Uhlenbeck random process. Then, obviously, it 
is quite impossible that recollision will not happen, since 
this would mean that symmetric Brownian path R(t) for 
the whole time t > remains confined below the straight 
line vt with negative slope v < 0. In opposite, the rec- 
ollision is inevitable. Consequently, its probability is not 
less than the probability that 2Vq < v < Vq < 0, that is 
(at nearly equilibrium situation) it has at least the order 
on ^Jm/M. 

Such estimates prompt that probabilities of multiple 
recollisions also are substantial. Due to all the recol- 
lisions, relatively slow atoms are accelerated by piston 
when it in its turn is accelerated by dominating relatively 
fast atoms. In other words, slow atoms "dress the pis- 
ton in a coat" thus randomly changing its effective mass. 
Seemingly, this must imply some fluctuations in relax- 
ation time and diffusivity of the piston. From the point 
of view of 0, @ this effect associates with that space- 
drift velocity of probability of some collisions chain is 
the center-of-mass velocity of particles taking part in this 
chain. Whether low-frequency fluctuations in diffusivity 
of Brownian particle wandering in one-dimensional ideal 
gas really exist is the question for separate consideration. 
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